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Abstract. This article examines the inverse problem for a lossy quantum 
graph that is internally excited and sensed. In particular, we supply an algo- 
rithmic methodology for deducing the topology and geometric structure of the 
underlying metric graph. Our algorithms rely on narrowband and visibility 
measurements, and are therefore of considerable value to urban remote sens- 
ing applications. In contrast to the traditional methods in quantum graphs, 
we employ ideas related to algebraic and differential topology directly to our 
problem. This neatly exposes and separates the impact of the graph topology 
and geometry. 



1. Introduction 

This article opens a new line of inquiry into the structure of waves on metric 
graphs. In particular, we are interested in algorithms for recovering the metric 
graph structure from measurements of propagating waves taken at unspecified lo- 
cations within the graph. We recover this structure in two stages: by first recovering 
the topology, and then the geometry. The algorithms for each stage draw upon a 
related circle of ideas in algebraic topology; in particular, we use homological tools 
to extract topology and cohomological tools to extract geometry. To validate the 
effectiveness of our approach, the algorithms have been implemented in software 
and tested on simulation data. |14j 

1.1. Historical development of the theory. The study of linear second-order 
differential operators on metric graphs has a long and venerable history, which is 
discussed nicely in the survey articles |26j [4] |15j . which have extensive references. 
The study of metric graphs with self-adjoint linear operators began with a short 
paper by Pauling |36j examining the spectral properties of aromatic compounds. 
This paper was followed by a number of others that refined Pauling's approach 
(for instance [40] , and of which [38] is a brief summary) . Because of the physical 
chemistry focus, a metric graph paired with a self-adjoint second-order linear op- 
erator became known as a quantum graph. Because the spectum of a self-adjoint 
operator is what is measured in quantum mechanics, the study of quantum graphs 
has been heavily focused on their spectral properties. The first mathematically 
rigorous treatment of the spectrum of a quantum graph was by Roth in [39 , which 
introduced the use of trace formulae. Trace formulae permit the spectrum of a 
quantum graph to be treated as a single algebraic object. 
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After Roth's initial work, an increasingly sophisticated literature on spectral 
properties grew up with trace formulae as the primary tool. As the spectral prop- 
erties themselves are not of direct interest here, we refer the reader to the excellent 
survey [15] and will call out a few interesting articles to highlight the historical 
development of the field. Spacing between eigenvalues in regular and equilateral 
graphs was treated by [20] and [34] . That the eigenvalues of a quantum graph are 
generically simple was shown by |12j , making use of a generalization of the theory of 
partial differential operators on manifolds. This simplifies some of the treatment of 
inverse spectral problems. Post |37) examines some general properties of the spec- 
trum of quantum graphs, and relates them to combinatorial Laplacians. Finally, 
Parzanchevski 35 showed that two different quantum graphs can have the same 
spectrum. 

In addition to being interesting in its own right, Parzanchevski's article 35i also 
demonstrates the kind of limitations present in recovering the underlying metric 
graph from its quantum graph spectrum. In other words, the inverse spectral prob- 
lem [8] cannot be solved completely. In line with spectral methods, which concern 
globally-valid solutions to linear equations on the quantum graphs, both inverse 
spectral problems and inverse scattering problems concern the extraction of met- 
ric graph structure from global measuments. The global nature of theses inverse 
problems is physically helpful, since quantum graphs have been used to study very 
small objects: it is difficult to excite a single atom in a molecule, usually one must 
settle for imposing global, external excitation. 

External excitations lead to scattering problems 31 , which in the case of quan- 
tum phenomena can lead to rather complicated dynamics for particles confined to a 
graph |24j . [231 . In particular, eigenfunctions of the Laplacian on a quantum graph 
can have widely separated peaks that are irregularly dispersed in the graph. This 
permits particles to tunnel in an apparently random fashion from one portion of 
the graph to another, and has therefore been a good model of quantum chaos. The 
examination of general scattering problems leads to interesting questions strictly 
outside of the quantum mechanical context. For instance, Flesia et at [11] ex- 
amined how the excitation of a regular graph with random speeds of propagation 
along it supports (or does not support) localized excitations. Indeed, they found 
that in classical wave propagation, the peaks of eigenfunctions tend to be more 
widely distributed than in the quantum case. 

Inverse spectral problems on quantum graphs became popular following the pa- 
per of Gutkin and Smilansky [17] . which was made somewhat more precise by 
Kurasov and his collaborators [27] , [29] . Briefly, these articles showed how to solve 
a generic class of inverse spectral problems when the linear operator is the Lapla- 
cian. The nongeneric problems clearly must have symmetries, as [5] explains in 
detail, and indeed Parzanchevski's examples of isospectral graphs are highly sym- 
metric. For other operators, there are partial results, of which [T] and [IS] are 
typical. For compact graphs in which the self-adjoint operator is the Laplacian, 
the spectrum determines the total length of the graph (sum of all edge lengths), 
the number of connected components, and the Euler characteristic (and hence the 
number of loops) [28] . 

Clearly, the topology of the underlying metric graph plays an important role in 
both the spectral structure and the solution to inverse problems. That connection is 
of particular interest in this article, as we show explicitly how algebraic topological 
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invariants of the graph impact the spectrum. For instance, in |13j . an index is 
computed that detects traveling wave solutions whose orbits are closed but not 
periodic when the self-adjoint operator on the graph is the Laplacian. (By "closed 
but not periodic," we mean traveling wave solutions that return to their initial 
starting point, but are going in the opposite direction at that starting point.) This 
index is topological, and depends only on the Euler characteristic of the graph. 

1.2. Main contributions of this article. In this article, we look at quantum 
graphs from a substantially different viewpoint, motivated by applications in re- 
mote sensing. In particular, while we are interested in inverse problems (finding 
the metric graph structure from measurements), we are interested in gathering lo- 
cal information only. To fix ideas, consider a dense urban environment, in which 
wave propagation can only occur in narrow channels (along the roads), and is 
obstructed by many buildings. Within this environment are placed an unknown 
number of transmitters and receivers at various locations. Given that the trans- 
mitters can be distinguished from one another by the receivers, how much of the 
geometric and topological structure of the network of propagation channels can 
be recovered? Since the propagation channels are narrow, it is sensible to assume 
that wave propagation happens only along the edges of a graph. The vertices of 
this graph represent junctions between propagation channels in the usual way. Of 
course, the resulting second-order operator is the wave operator, which leads to a 
quantum graph formulation of this problem. (This "thin-wire" approximation has 
been examined asymptotically in [35], [25], [55], and [35].) Unlike the inverse spec- 
tral and scattering problems, the measurements are from internal locations, rather 
than at a distance and concern local features of the data. Indeed, the continuity of 
the eigenfunctions is what will permit us to infer global features from these local 
measurements. 

There are three main reasons why the traditional quantum graph analytic tech- 
niques would be impractical in a remote sensing application: 

(1) In order to collect enough spectral data to exploit any of the inverse prob- 
lem solutions found thus far, extremely wide-band receivers and transmit- 
ters would be needed. In the case of urban remote sensing, the resulting 
frequencies are usually not available for imaging purposes. 

(2) Propagation losses can be substantial in remote sensing, which usually re- 
sults in a strong localization of the signal to the vicinity of each transmitter. 

(3) The computational complexity of the algorithms resulting from the inverse 
spectral problems is unknown, and may hamper practical application. 

These difficulties suggest that a better approach is to exploit transmitter visibility 
information (ie. tabulate which transmitters are visible by which receivers). As will 
be shown, this permits the topology of the underlying graph to be computed algo- 
rithmically and efficiently. We exhibit a novel algorithm for preprocessing visibility 
information so that it is suitable for computing the underlying topology. Coherent, 
local measurements of the waves permit the geometry to be recovered, under sim- 
ilar genericity conditions to |17j . Additionally, we obtain a new characterization 
theorem that relates the space of eigenfunctions to the geometric and topological 
structure of the graph. 

The tools we employ are those in the emerging field of applied algebraic topology. 
In particular, we will employ the Nerve Lemma [6] to determine the topology from 
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visibility information, the Whitney embedding theorem for stratified spaces [33j . 

and some ideas from the cohomology of constructible sheaves. 

It should be noted that the most similar approach that we have found in the 
literature is that of the recent work by Caudrelier and Ragoucy [9], which uses 
reflection-transmission algebras to perform computations on internally-sampled sig- 
nals. We recover their results using a sheaf-theoretic framework, which cements the 
connection of their method to topological invariants. 

1.3. Outline of the article. We begin in Section[2]by giving a precise statement 
of the problem we aim to solve, as well as stating the key results that we prove in 
later sections. In particular, we observe that there exists (Theorem[3]) a finite cover 
by regions coming from thresholded signal levels. This is an enabling result for 
an algorithm (Algorithm [9]) for extracting the graph topology from measurements. 
We explain this algorithm in Section G3 In Section |4j we demonstrate that once 
the topology has been found, narrowband signals suffice to compute the geometry 
(Theorem 1361 which has an inductive, algorithmic proof). In support of this result, 
we prove a characterization (Theorem 1231) that shows how the space of solutions 
to quantum graphs is completely, and explicitly determined by the geometry and 
topology of the graph. Finally in Section [5] we discuss the results. 



2.1. Review of metric graph definitions. The underlying metric graph struc- 
ture that will be used in this article is essentially the one initially defined in (and 
later examined in [3], [2]), but to support the examination of scattering problems 
we relax the compactness hypothesis. 

Definition 1. A metric graph is a metric space X such that each x £ X has a 
neighborhood U x isometric to a degree n x star-shaped set: 



We can associate a combinatorial structure of edges and vertices (nonuniquely) 
to a metric graph. Precisely, if V C X is a discrete subset such that V contains 
all points of X whose degree is not 2, we call V(X) a set of vertices. The set of 
path-connected components of the complement of V (each of which is isometric to 
an open interval) is called the set of edges E(X). Those elements of V(X) that lie 
in the topological closure of an edge e are called the endpoints of e. (There may 
be only one endpoint in any given edge.) (We remark that the choice of vertex 
and edge sets will in general effect the quantum graph structures we obtain, but we 
choose vertex conditions for which this nonuniqueness is immaterial.) 

We assume that it is always possible to find a finite set of vertices and edges. 
We call X in this situation a finite metric graph. With this structure, an edge that 
contains a degree 1 vertex will be called a closed edge. An edge whose closure in X 
is not a circle, has exactly one endpoint, and has infinite length will be called an 
open edge. 

Of course, this assumption replaces the compactness assumption in [46] and 
others. 

The primary contribution of [3] was the formulation of a Laplacian operator 
suitable for finite metric graphs. For a piecewise smooth function / : X — > C, this 
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operator is the bounded, signed measure given by 

(1) A/(x) = 0(x)+^a x (/)5 I , 

where cr x (f) is the sum of (outgoing) directional derivatives of / at xQ Notice in 
particular, that the second derivative operator appearing in the first term of (fT]) 
is insensitive to the parametrization at x on the interior of an edge. Secondly, the 
sum of directional derivatives vanishes everywhere except at finitely many points 
since / is piecewise smooth. Since the Laplacian is a measure, the fact that the 
second derivatives are not defined at the vertices presents no difficulty 

2.2. Lossy quantum graphs and their fundamental solutions. Since we are 
interested in narrowband, lossy wave propagation, we will study the fundamental 
solution to the Helmholtz equation want to examine solutions to 

(2) Au + k 2 u = S y , 

where the wavenumber k may be complex (to indicate lossy propagation), and 
y G X is the transmitter location. We will be interested in solutions that satisfy 
the "Kirchoff" conditions |26) : 

• u is continuous on X and 

• at each point x £ X, the sum of the derivatives of u, o~ x (u) = 0. 

Definition 2. The pairing of @ and the Kirchoff conditions will be called a lossy 
quantum graph. 

Clearly, along the interior of an edge, one has the general solution to <j2j) which 
is the superposition of two traveling waves: 

(3) u(x) = Cl e lkx + c 2 e- lkx . 

The non-unique choice of vertex and edge sets does not play a significant role in 
the solutions of a lossy quantum graph when Kirchoff conditions are used. Suppose 
u is the solution to a lossy quantum graph. If we consider another quantum graph 
structure in which there is an additional vertex within an edge, u will automatically 
be a solution of this new structure. Likewise, any solution of this new structure 
will be a solution of the original structure. This latter fact is most easily seen as a 
consequence of edge collapse (Lemma [501) • 

A fundamental solution (f2]) of a quantum graph X can be obtained as a super- 
position of solutions of the lossy quantum graph on X — {y}. It can be interpreted 
as an external source scattering problem for the graph with y removed. 

2.3. Thresholding and covers of contractible regions. In a lossy quantum 
graph, the energy from a signal source tends to be concentrated around its loca- 
tion. From an urban signal propagation perspective, this means that only receivers 
nearby a given transmitter will be able to decode its transmissions correctly. This 
suggests that visibility within a lossy quantum graph could be useful for determin- 
ing which receiver locations are near which transmitters. System designers typically 
use signal level thresholds to determine whether decoding can proceed successfully, 
so we employ a thresholding approach to construct visibility regions for a given 
signal source's location. Precisely, the region where the signal level is above a given 



Caution: our sign convention differs from [3]. 



6 



MICHAEL ROBINSON 



U 



2 








X 



L 



Figure 1. Amplitude of interfering traveling waves 



threshold will be a contractible subset of X when the threshold and signal loss are 
large enough. 

The central result is that there exists a collection of locations of signal sources 
that can be thresholded appropriately to give a collection of contractible visibility 
regions that cover the graph. 

Theorem 3. Suppose X is a lossy quantum graph with wavenumber k = k' + ia 
for k',a > 0. (a is the loss coefficient, with larger values corresponding to higher 



(1) For each y G X , there is a choice of a (say a y ), and a threshold T y > such 
that for the fundamental solution u y to @, the visibility region U y (a y ) = 
{x £ X\\u y (x)\ 2 > T 2 } is contractible. 

(2) The collection {U y (a y )} forms a cover of X . 

(3) Suppose that all vertices of X are contained in the compact set K . There is 
therefore a finite subcollection {yi, ■■■,y n } such that {U Vl (a m ), U yn (a yn )} 
covers X n K. Let a — max{Q Sl ,...,a 9n }. Then the resulting collection 
{U Vl (a), U Vn (a)} covers X n K and consists of contractible sets. 

In order to prove Theorem [3l we need to study the interference of two lossy 
traveling waves on a line segment. The requisite information is captured in the 
following Lemma. 

Lemma 4. Consider < x < L and two traveling waves u(x) — e lkx + ce~ lk ( x ~ L ^ 
with k = k' + ia as in the statement of Theorem^ Suppose that c = Ye~ aL , where 
r is the reflection coefficient at x — L. Note that the amplitude of the right- going 
wave at x = L is e~ aL . For any fixed T, k' , and L, there exists an a > and a 
T > such that {x € [0, L]\\u(x)\ > T} is connected and contains zero. 

Note that T is strictly positive. 

Proof. We proceed by computation... 

u = e - ax e lk ' x + ce ^-L) e -ik{x-L) 



e - ax cosk'x + ce a{x - L) cosk'(x -L))+i (e- ax sink'x - ce a(x - L) sink'(x 



loss.) 
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Thus, after a small amount of manipulation, 

(4) M 2 = e- 2ax + c 2 e 2a{x ~ L) + 2ce~ aL cos {2k 1 x - k'L) . 

Notice in particular that the fast spatial variation is captured by the cosine term. 
The location of the "first sidelobe" (see Figure [1} of this expression is the left most 
location of the maxima of the cosine, namely 

717T L 

XFSL = 2k' + r 

where n is the smallest integer such that Xfsl is positive. 
An upper bound for \u\ 2 is its envelope, namely 

(5) M 2 < e~ 2ax + C 2 e 2a{x - L) + 2ce~ aL . 
At the first sidelobe location, the value of | tt | 2 is 

\u{x FSL )\ 2 = e- aL (e-T* 21 +c 2 e^ + 2c) 

= e~ aL (e - ™ +T 2 e £ ^- aL + 2T e - aL ) , 

which can be made less than the value at 

|w(0)| 2 = 1 + T 2 e- 3aL + 2Te- 2aL cos(-fc'Z) w 1. 

by taking a large enough. Therefore, we'd like to take T = \u(xfsl)\- 

Now the upper bound (O is symmetric about x — L ~ l ° sc ( so the only thing that 
could spoil the connectedness of the super- level set of T is large values of \u\ near 
L. Hence, we want 

L < the last sidelobe location before the minimal envelope 

nir L . 
< L — log c — — — — — (an underestimate) 

L vm 
2 < " l0gC ^2F- 

L(±-a) < -logr-^, 
\2 J 8 2k' ' 

which is clearly satisfied for large enough a. □ 

Remark 5. The wavenumber k' is fairly high for practical imaging systems (roughly 
tens of waves per meter), so one typically spatially averages the measurements to 
avoid aliasing. This means that instead of thresholding the amplitude ((4|) 7 we thresh- 
old the envelope ([5]). This results in a tolerance of lower loss, and also results in 
larger visibility regions. 

Now, we can present the proof of Theorem |3j 

Proof. (1) This is immediate from repeated application of Lemma |4j and the 
conversion of each edge incident to a source y into an externally excitated 
edge. The source is at y, and the reflection coefficient T captures effects 
from the rest of the graph as well. Notice that more loss will make T 
smaller, and thus a threshold that satisfies Lemma |4] will continue to work. 
(2) That the visibility regions form a cover is immediate from the fact that 
each region obtained in Lemma 0] contains zero. 



s 



MICHAEL ROBINSON 



(3) The important point is that for larger a, the sets U y (a) can be increased in 
size, or even maintained at the same size (on a single edge) by the selection 
of a threshold. In any event, it is apparent that the since the degree of any 
vertex is finite, a single threshold per source is easily found that results in 
cover by contractible sets even when the loss is increased. 

□ 

Remark 6. Theorem^ is an existence result only: many configurations of trans- 
mitters will fail to provide a cover upon thresholding. However, by adding finitely 
many more transmitters to areas of low transmitter density, a cover by contractible 
sets can be obtained. 

3. Topology computation 

In this section, we make use of the visibility regions of the previous section 
(Theorem [3]) to deduce the topology of a quantum graph. Our primary tool is the 
Nerve Lemma, generally attributed to Leray [6]. If one has a cover of the graph 
by contractible sets (as in Theorem [3]), for which all finite intersections between 
elements of the cover are also contractible, then one can construct a simplicial 
space that is homotopy equivalent to the original graph. This nerve construction is 
easy to implement, and for which there are computational topology tools available, 
for instance |21) . 

Unlike the work of previous authors, the algorithm presented here (in Section 
13. 2p only relies on coarse discriminations, rather than on a detailed examination 
of the spectrum. In particular, the only assessment that is needed is whether 
the intersection between a number of visibility regions is connected or not. It is 
not immediately clear that this determination is possible. However, according to a 
generalization of Whitney's celebrated embedding theorem [30] , there is a topology- 
preserving embedding of the graph into an appropriately high-dimensional signal 
space, which we construct. The result is that if the intersection of several visibility 
regions is disconnected, then its image is also disconnected in the signal space. 

The main idea of the proof of the algorithm's performance centers around the 
Mayer- Vietoris sequence for homology, which is a well-understood tool in algebraic 
topology. We refer the reader to the excellent treatment in [18] for details. 

3.1. Uniqueness of signal response and the detection of connected com- 
ponents. If the visibility regions found in Theorem [3] formed a good cover, we 
would be able to use the Nerve Lemma to recover the homotopy type of the un- 
derlying metric graph. However, the visibility regions will usually not form a good 
cover, often in the way shown in Figure [2] However, in a graph, the intersection 
of two contractible sets can fail to be contractible only by being disconnected. In 
this section, we show how to detect disconnectedness of sets by using three or more 
independent solutions to the quantum graph. 

If there are m signal sources, the complete set of measurements we obtain at a 
receiver located at x is the value of the piecewise smooth function 

P(x) = (u 1 (x),u 2 (x),...,u m (x)). 

For brevity, we say that P is the signal profile, taking locations in X to signal space 
S. 
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If the signal profile is topology-preserving (an embedding), then U C X is dis- 
connected if and only if its image under P is disconnected. Better, we could just 
look at the image (which are measurements taken from every point) and deduce 
the homeomorphism type of the graph. (Indeed, homeomorphism type is stronger 
than homotopy type.) However, this is unrealistic: normally we have to sample 
the receiver locations. This is the reason for using a coarser construction based on 
refining the visibility regions into a good cover. 

We can obtain the following result (appearing for manifolds in 16 ) using a 
stratified transversality result in [U] . 

Theorem 7. If the dimension of the signal space m (the number of transmitters) 
is greater than 2, the signal profile is generically (almost always) an embedding. 

Proof. First, recall that an embedding consists of an injective, piecewise smooth 
function whose derivative (where defined) is nonvanishing, by the inverse function 
theorem. By direct computation along the edges, the derivative of a nontrivial 
solution to the Helmholtz equation is nonvanishing. 

It is also therefore the case that a metric graph is an (a)-regular Whitney strat- 
ified space, ((a)-regularity is essentially trivial in this case since the only strata are 
vertices and edges, and vertices have trivial tangent space. See [44] for a precise 
definition.) 

Rather than working with P directly, we work with an augmented version 
F : C 2 (X) -> C 2 (X) xSxS, 

given by 

F{x,y) = (x,y,P(x),P(y)), 

where C 2 (X) = X x X — Ax is the space of all nonoverlapping pairs of points on 
the graph. (Ax is the diagonal in X and consists of pairs of points at the same 
location.) We note that C 2 is an (a)-regular stratified space, being the product of 
two such spaces (removing the diagonal causes no difficulty). Since F is evidently 
an injective immersion, the image of F is also (a)-regular. (We note that without 
the factor of C 2 in the codomain, F fails to be injective whenever P is not injective; 
this disrupts the stratification of the image.) 

P fails to be injective if and only if the image of F intersects D = C 2 x As. 
We use transversality to determine whether that intersection is empty or not. The 
transversality result of [H] states that F is generically transverse to D. (We actually 
perturb the submanifold D, not the image of F. This amounts to perturbing not 
the recieved signals, but our ability to determine if the signals from two points 
agree. We are looking for an empty intersection between the stratified set and the 
diagonal. What we get then is stability of this empty intersection.) 

When F is transverse to D, then the intersection between any stratum of the 
image of F and D has dimension no greater than 

dim C 2 + dim D - dim (C 2 x S x S) = 2dim X + (2dim X + m)- (2dim X + 2m) 

— 2dim X — m = 2 — m. 

When this bound on dimension is negative, there is no intersection under generic 
perturbations. Hence, if the number of independent signals measured is greater than 
2, the signal profile is generically injective, hence generically an embedding. □ 
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Obstacle 



Figure 2. Intersection of two contractible coverage regions may 
be disconnected. (According to Theorem another transmitter 
would be needed to disambigute the components of the intersec- 
tion.) 

3.2. Coverage refinement. The visibility regions associated to each transmitter 
can be made contractible by appropriate thresholding. However, this usually is not 
enough to make the resulting cover into a "good" cover, and thereby recover the 
topology. In a graph, the only obstruction is that intersections between coverage 
regions may be disconnected. In Section [3. 1[ we showed that disconnected compo- 
nents could be detected and identified from the received signals. In this section we 
exhibit an algorithm that exploits connected component discrimination to refine 
the coverage regions into a good cover. 

In this section, we use homology with a fixed field F as the set of coefficients. We 
will therefore suppress the coefficients from the notation. This has the advantage 
that in graphs, we can identify acyclic sets with contractible ones [21] . 

Theorem 8. Suppose that U = {Ui, U%, ...t/ n } are open sets forming a good cover 
of a subgraph of X and that W is a contractible subset of X for which 

• W <£l)U, and 

• there are U^, Ui 2 , ... € U such that W H Ui 1 n Ui 2 n ... is disconnected. 
Let V = W n (UK), then 

(1) each path component of V is acyclic, and 

(2) there exist open neighborhoods XJ[ of ' f/j — V for each Ui € U and an open 
neighborhood W' ofW — V that satisfy 

(a) U{{7 4 '} and W are disjoint open sets, 

(b) each XJ[ and W consists of acyclic path components , and 

(c) the collection of the path components of V , W , and the U[ forms a 
good cover. 

It is straightforward to make use of Theorem [8] to devise an iterative algorithm 
for refining a set of contractible transmitter coverage regions V = {Vx, V%, V n }. 
We have implemented this algorithm in software and tested it against computer 
simulations. [14] 
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Algorithm 9. We define good covers Vi, V2, V n inductively: 

• Base step: V\ — {V±} is a good cover by assumption. 

• Induction step: Assuming Vk is a good cover, we apply Theorem [3| with 
W = Vk+\ and U = Vk- We then define Vk+i to be the good cover that we 
obtain as (2)(c) in Theorem^ 

We begin by proving a technical lemma that we primarily use to relate the 
homology of a graph to a subgraph. This is an easy application of the Mayer- 
Vietoris sequence. 

Lemma 10. Suppose that f : X — > Y is a continuous injection from one graph 
into another and that Y — /(A) is a disjoint union of open, acyclic sets. Then f 
induces an injection H\(X) — > Hi(Y). 

Proof. Simplify this proof!! First, observe that there is an open neighborhood 
U of f(X) that is homotopy equivalent to f(X). This neighborhood consists of the 
union of f(X) and some subintervals of Y — f(X) (taken from the ends of the edges 
in Y — f(X) that are adjacent to f(X)). Then we may consider the Mayer- Vietoris 
sequence 

-> ffi(some intervals) -> Hi(U) ® H X {Y - f{X)) -> H X {Y) -> ... 

since U U int (Y — f(X)) = Y. By assumption U has the same homology as f{X) 
and X, so this exact sequence can be written 

-» ffi(some vertices) -> Hi(X) - > H X {Y) -> ... 
A discrete space always has trivial Hi, whence /* must be injective on H%. □ 

The proof of Lemma [8] relies on the topological dimension of X being not greater 
than 1. Indeed, the Lemma is false in dimensions 2 and higher. 
We now proceed with the proof of Theorem [8l 

Proof. Begin by considering the Mayer- Vietoris sequence 

-> Hx(V) -> ffi(UW) © Hi(W) — Hi(W U (UK)) H (V) -4 ... 

Note that H\(W) is trivial by assumption, and by Lemma fTOl the map * above is 
injective. Hence each component of V is acyclic, establishing (1). 

Define B t — Ui Pi (V — V) = Ui f) dV, which consists of the boundary points of 
V that lie in fy. Similarly, define C = WC\dV. 

Observe that C is disjoint from all of the Bf. 

cniuB,) = (w n dV) n (u([/j n dV)) 
= wn(uUi)ndV 
= v n dv = 



since V is open. 

Further, dV is finite, since A is a metric graph. 
Define 

(6) S = min d(x, y). 

x^y 
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Since C U (U{B,}) is finite and contains more than one point, 6 > 0. However, an 
upper bound for S is 

(7) 6 < m£{d(x, y)\x e W - V, ye (UU) - V}. 

Define neighborhoods 

W' = (W-V)U I (J B s/3 (x)) , 

\x£C J 

and for each i, 

Ul = (U i -V)ul\J B s/3 (x)\ , 
\xeB z ) 

where B r (x) is the open ball of radius r centered at x. Evidently by (0, W'nU! = 
for all i, which establishes (2)(a). By ©, 

W' ~ W - V, and U-c^Ui- V. 

Notice that H\{W — V) and H\(Ui — V) are trivial by Lemma [TU1 consider the 
inclusion of W — V — > VF. We just proved that V is a disjoint union of acyclic 
components, and it is evidently open. This inclusion induces an injection on H±. 
However, W is acyclic, so this implies that H±(W — V) must be trivial. The same 
reasoning works for U% — V, which establishes (2)(b). 

To show (2)(c), it remains to establish the following statements, which are im- 
mediate from the construction of W' and U-: 

• Components of V and W have acyclic intersections since they're precisely 
half-open intervals, 

• Components of V and U( have acyclic intersections by the same logic, and 

• The collection of components of {[//} form a good cover of (UU) — V. 
Observe that disjoint intervals were attached to each {/, in the construction 
of XJ[ such that if an interval 1{ was attached to Ui and Ij was attached to 
Uj, with Ii [~l Ij ^ 0, then I t = I r 

□ 



4. Geometry computation 

We now address the problem of obtaining the geometric structure of a quantum 
graph by using a few narrowband measurements taken at points internal to the 
graph. The central result of this section shows how a single solution over the entire 
graph determines the geometry up to phase ambiguity. It should be emphasized 
that this result is obtained under the assumption that the topological structure 
of the underlying metric graph is known, either as a homeomorphism type (which 
is preferable), or more likely as a nerve. In either case, we assume that we are 
in possession of a 1-dimensional simplicial complex that describes the underlying 
topological space structure of the quantum graph. 

Since the primary effect of geometry is on the phase of signals, the measurement 
of geometric information in a quantum graph requires coherent, time-sensitive pro- 
cessing. From a signal processing perspective, the resulting signal-to-noise require- 
ments are considerably more demanding than those required to obtain visibility 
information. As a result, the measurements of geometry in a lossy quantum graph 
will tend to be useful over a fairly small region. We therefore compute geometry 
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from local measurements, which are the least contaminated by loss. Additionally, 
we want to operate passively, without reference to known transmitter locations. 
Therefore, without substantial loss of generality, we consider homogeneous solu- 
tions rather than fundamental solutions. 

As an aside, the calculations given in this section are actually still valid for 
homogeneous solutions of lossy quantum graphs. These solutions correspond to 
waves incident from outside the graph. (See the discussion following Question [24] 
for details.) 

The space of solutions to a given quantum graph is finite dimensional, and de- 
pends on both topology and the geometry. In Section 14.31 we compute this di- 
mension from the graph structure. In more traditional language, we obtain the 
dimension and structure of the eigenspaces of the graph Laplacian operator. (There 
is a delicate interplay between loss and the dimension of these eigenspaces: signal 
loss eliminates resonance phenomena that are extremely useful in inverse spectral 
methods for detecting geometric and topological features.) Sheaves are the appro- 
priate mathematical tool for connecting local information to global behavior, so we 
develop sheaf theoretic computational tools for quantum graphs. We obtain sheaf- 
theoretic proofs of certain graph operations that preserve quantum graph structure, 
which allows great computational simplification. 

We give a brief introduction to the key ideas of sheaf theory in Section 14.11 For 
a more detailed exposition, we refer the reader to Appendix 7 of [19] and to [7]. 
Building on this introductory material, we give an explicit definition of the sheaf 
structure to be used on a quantum graph in Section 14.21 The global structure of 
this sheaf is computed in Section 14.31 (Theorem [23)) , which explicitly demonstrates 
the topological and geometric dependence of the space of solutions. As an aside, 
we note how this recovers an inverse spectral result of previous authors. Finally, 
in Section 14. 4| we show how combining the topology and a single solution provides 
detailed geometric information about the graph. 

4.1. A Brief introduction to sheaf theory. A sheaf is a mathematical tool for 
storing local information over a domain. It assigns some algebraic object, a vector 
space in our case, to each open set, subject to certain compatibility conditions. 
These compatibility conditions are of two kinds: (1) those that pertain to restricting 
the information from a larger to a smaller open set, and (2) those that pertain to 
assembling information on small open sets into information on larger ones. What 
is of particular interest is the relationship of the global information, which is valid 
over the entire graph, to the topology of that graph. Additionally, sheaf theory 
identifies classes of transformations of the underlying graph that preserve the global 
information. These transformations permit us to simplify the graphs with no loss 
of generality. 

4.1.1. Elementary definitions for sheaves. In this section, we follow the introduction 
to sheaves given in Appendix 7 of |19) , largely for its direct treatment of sheaves over 
tame spaces. For more a more general, and more traditional approach, compare 
our discussion with [7] . 

Definition 11. A presheaf F is the assignment of a vector space F(U) to each open 
set U and the assignment of a linear map p(j : F(U) — > f or each inclusion 

V C U. We call the map p(j the restriction map from U to V. Elements of F(U) 
are called sections of F defined over U. 



14 



MICHAEL ROBINSON 



Definition 12. A sheaf T is a presheaf F that satisfies the gluing axioms: 

• (Monopresheaf ) Suppose that u e F{U) and that {U\, U2, ■■■} is an open 
cover of U . If p\jU — for each i, then u = in J-{U). Simply: sections 
that agree everywhere locally also agree globally. 

• (Conjunctivity) Suppose u G F(U) and v G F{V) are sections such that 
Pu nv u = py nV v. Then there exists awe T{U U V) such that p^ uV w = u 
and pJjyjyW — v. In other words, sections that agree on the intersection of 
their domains can be "glued together" into a section that is defined over the 
union of their domains of definition. 

Example 13. Standard examples of sheaves are 

• The collection of continous real-valued functions C {X, R) over a topological 
space X. In this case, the sections defined over an open set U is {/ : U — >■ 
IR|/ is continuous}. 

• The collection of locally constant functions, which essentially assigns a con- 
stant to each connected component of each open set. 

In contrast, the collection of constant functions does not form a sheaf. Suppose 
u and v are distinct constants defined over disjoint sets U and V. The sheaf axioms 
would indicate that since their domains of definition (U and V ) are disjoint, then 
there should exist a constant function defined over Ul)V that restricts to each. This 
is of course impossible, since such a function is only locally constant. 

We now turn to the problem of understanding the effects of graph operations 
on sheaves. There are six famous operations on sheaves that are important in the 
general theory, but only two of them (cohomology and direct images) play a role in 
this article. 

4.1.2. Cohomology. We can recast the conjunctivity axiom as testing if Pij nV u — 
Py nV v is zero or not, rather than checking for equality. This can be viewed as 
looking at kernel of the linear map d : J r (C/)©J r (V^) — > F(UC\V) given by d(x,y) — 
Pij nV x — py nV y. Indeed, all of the elements of the kernel of such a linear map 
correspond to the agreement of sections on U n V. 

On the other hand, the monopresheaf axiom indicates that the preimage of zero 
under the map d corresponds to the restriction of these glued sections onto each of 
U and V. Indeed, any nonzero element of the image of d cannot be a section over 
UUV. 

These two points motivate a computational framework for working with sheaves, 
called the Cech construction. 

Definition 14. Suppose J 7 is a sheaf on X , and that U = {U\, U2, ■■■} is a cover of 
X. We define the Cech cochain spaces C k (U;J-) to be the direct sum of the spaces 
of sections over each k-wise intersection of elements in U. That is 

C k {U;F) = ($F{U ll C\U l2 C\...C\U lk ). 

We define a sequence of linear maps 

d k : C k {U;F) -> C k+1 {U;F) 

by 

h-\- 1 

d k (a)(U 1 , U 2 , U k+1 ) = ^(-l) l pTn:'^ k T +la ^ n -' <- n ^+1). 

i=0 



INVERSE PROBLEMS IN GEOMETRIC GRAPHS USING INTERNAL MEASUREMENTS 15 



where the hat means that an element is omitted from the list. Note that these fit 
together into a sequence, called the Cech cochain complex: 

0^C Q (U;T) ^{U-T) — ... 

A standard computation shows that dk ° dk-i — 0, so that we can define the k-th 
Cech cohomology space 

H k (Li;J 7 ) = ker dk/ image dk-i- 

The H k apparently depend on the choice of cover U, but for good covers (much 
as in the Nerve Lemma), this dependence vanishes. Leray's theorem for sheaves 
states that H k (U\ J 7 ) is the same for each good cover. Indeed, it then depends only 
on X. So we write H k (X;T) — H k (U;F), the sheaf's cohomology in the case that 
U is a good cover. 

A little thought about good covers on graphs reveals two important facts: 

• for a metric graph X, H k (X; F) = for k > 1, and 

• H°(X; J 7 ) is isomorphic to the space of global sections J-(X). 

The first fact follows immediately from a covering dimension argument. The 
latter fact comes from our construction, and that the image of <i _1 (not d inverse!) 
is zero. This suggests a computational way to examine solutions to quantum graphs 
by way of sheaf theory: we construct a sheaf that has local sections being the locally 
valid solutions to the quantum graph, and then compute its cohomology. What is 
especially valuable about this approach is that there is additional information in 
J 7 ). This additional information plays a useful role in Section POl and detects 
resonance phenomena in quantum graphs. 

By analogy with the Mayer- Vietoris sequence for homology, there is a Mayer- 
Vietoris sequence for sheaf cohomology. It is given by the following theorem (also 
in [7] in a variety of forms) : 

Theorem 15. Suppose that A, B are two open subspaces of a graph X that cover 
X, and that J- is a sheaf over X. Then the following Mayer- Vietoris sequence is 
an exact sequence: 

...^H k {X;T) — ^— >• H k (A; T) © H k (B\ T) — ^ H k (A n B; T) — S —> H k+1 (X; T) — > ... 

In this sequence, r comes from restriction maps in the obvious way, d is the com- 
position of restriction maps and a difference: d(x, y) = p^ nB x — p^ nB y, and S is 
the connecting homomorphism. Note that notation has been abused above slightly: 
by H k (A;J r ) we mean the k-th cohomology of the sheaf J- restricted to subsets lying 
in A. 

4.1.3. Direct images. 

Definition 16. Suppose f : X — > Y is a continuous function between topological 
spaces. If J- is a sheaf over X , then the direct image of T through f is written f*T 
and is given by its action 

{UT){u) = Hr 1 (u)) 

on open sets U CF . Clearly, f*T is a sheaf over Y . 

The key question is when the cohomologies of T and f*J- agree. If they do 
agree, and the combinatorial structure of Y is much simpler than X, then it is 
usually much easier to compute the cohomology of f^T rather than of T . We 
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make extensive use of this in the sequel. This essential question is answered by the 
following result: 

Theorem 17. (The Vietoris mapping theorem, Theorem 11.1 in Chapter II of 
Suppose that f : X — > Y is a closed, continuous function between metric graphs, 
and T is a sheaf over X. If H k (f~ 1 ({y});J 7 ) for k > and each y £ 7 , then T 
and f*J- have isomorphic cohomology. 

4.2. Definition of a sheaf structure for quantum graphs. The collection of 
locally-valid, continuous solutions to a differential equation on R forms a sheaf: 
restricting the domain of validity of a local solution produces a new local solution, 
and local solutions can be joined if they agree on the intersection of their domains 
of validity. This construction also works on a quantum graph, resulting in the sheaf 
of excitations. 

Since solutions on a single edge consist of two propagating waves, we obtain 
a more explicit definition for the sheaf of excitations. We construct the sheaf of 
excitations as a direct image of a sheaf over a directed graph with two oppositely- 
oriented edges for each undirected edge of our original quantum graph. [55] Our 
construction permits a slight generalization over the sheaf of excitations in that the 
values of the solution may take values in any field rather than the complex numbers. 
This simplifies the notation, so we perform all further computations at this level of 
generality. The resulting general object is called a transmission-line sheaf. 

Proposition 18. Local solutions of a quantum graph form a sheaf, called the ex- 
citation sheaf. 

Sheaves provide a convenient way to collate solutions to differential equations. 
For instance, [43] discusses the situation in great detail for systems of linear partial 
differential equations. We sketch a proof here for our specific instance. 

Proof. The appropriate sheaf is a correspondence between open sets (of the topo- 
logical space structure of the graph) and the space of continuous functions over 
them, in which inclusion of open sets corresponds to restriction of functions. Solu- 
tions on an open set in the interior of an edge remain solutions when restricted to 
a smaller domain. Additionally, if solutions on two open sets within the interior of 
an edge agree (pointwise) on the intersection of these two open sets, then there is 
a solution defined on the union of these open sets. 

The vertex conditions require continuity of solutions and some additional con- 
straints. Suppose that s is a solution defined on U, an open set containing a vertex 
v. Clearly, s satisfies the vertex conditions at v. Now if V C U is another open 
set, restricting s to V will result in another solution, fn particular, if v £ V, then 
the s\V will still satisfy the vertex condition at v. Now if si and S2 are solutions 
defined on U\ and U2, and v lies in U\ n U2, then both Si and S2 will satisfy the 
vertex conditions at v. If si\(Ui D U2) = S2I (C^i H U2), then we can define a new 
solution s pointwise by 



This new solution clearly satisfies the vertex condition at v, and is otherwise still 




si (x) if x 6 U\ 
s 2 {x) if x G U2 



at solution. 



□ 
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The excitation sheaf is a complicated object, one in which the combinatorial 
structure is obscured. In particular, it is unclear how to compute its sheaf coho- 
mology, which encapsulates all of its associated global structure. Therefore, we 
construct a new sheaf called a flow sheaf with simpler structure that has isomor- 
phic cohomology. The decisive feature of this flow sheaf is that it is constructible 
with respect to the graph structure. This means that over edges or vertices, the 
sheaf structure is constant. 41 Only where edges attach to vertices does the sheaf 
structure change, and then only in a constrained fashion. 

A flow sheaf, however, is defined over a directed graph. We'll use a direct image 
construction to write the appropriate sheaf definition for quantum graphs, which 
are undirected. 

Definition 19. Suppose that ¥ is a field and that X is a directed graph. Given that 
X has the usual topology, let 1A = {U a , Vp} be a base for the topology of X where 
each U a is connected and contains exactly on vertex and each Vp is contained in the 
interior of a single edge. An W-flow sheaf on X is the sheafification of the following 
presheaf F , defined on U: 

• F(U a ) is a direct sum of copies F, one for each incoming edge into the 
unique vertex contained in U a , 



• ifVpdU a and Vp is contained in the n-th incoming edge, the the restriction 
map F(U a ) —> F(Vp) is projection onto the n-th copy of¥, 

• if Vp C U a and Vp is contained in the n-th outgoing edge, then there is 
a fixed ¥ -linear map F(U a ) — > F(Vp) depending only on the vertex v con- 
tained in U a and n (the outgoing edge). This collection of maps, one for 
each outgoing edge, is called the local coding at v and is denoted by <p n (v). 

The following proposition is immediate from the definition of a constructible 
sheaf. 

Proposition 20. Flow sheaves are constructible with respect to the skeleta of X . 

Now let us consider the more specific case of a direct image of a flow sheaf that 
is isomorphic to the excitation sheaf defined earlier. Suppose £k is an excitation 
sheaf with wavcnumbcr k over a quantum graph X . Let V and E denote the vertex 
and edge sets of X, respectively. Denote by L e the length of edge e G E. 

Construct a new graph Y (that is a directed metric graph), by replacing each 
undirected edge Vi <-> v-i of X by a pair of opposing edges v\ — > V2, v% — >■ v\. 
Let / : Y — > X be the map that sends each such pair of directed edges to the 
corresponding undirected edge of X, as shown in Figure [3l We construct a C-flow 
sheaf T over Y with coding maps tailored to work like the Kirchoff conditions, as 
follows. Consider a vertex v € V, with degree n m X. Label the incoming edges 
(in Y) {at,a%, ...a n } and the outgoing edges {61,62, — b n } so that /(a,) = f(h) for 
each i. Define the coding map by 



where zi is the value of a section of T restricted to the i-th incoming edge 

We call Gk — f*J~ a transmission line sheaf. Any sheaf constructed over an 
undirected graph as the direct image of an F-flow sheaf with coding maps satifying 



• F(Vp) = ¥, 



(8) 
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X Y 

Figure 3. Construction of a transmission line sheaf on an undi- 
rected graph from a flow sheaf on a directed graph 




Figure 4. Example of a transmission line sheaf 



(JSJ (but with F endomorphisms specified on the edges, rather than phases) will also 
be called an F-transmission line sheaf. 

Suppose J 7 is a transmission line sheaf on a graph X. A resonant edge is one 
in which the edge endomorphisms are identity maps. If the edge endomorphisms 
differ from the identity, we call the edge nonresonant. 

Proposition 21. The excitation sheaf Zk and the transmission line sheaf Gk ore 
isomorphic as sheaves. As a result, we may study the solutions to a quantum graph 
by computing the cohomology of the transmission line sheaf instead. 

Proof. One need only observe that the Kirchoff conditions are equivalent to the 
coding map §8§ employed at the vertices of X, by a straightforward algebraic trans- 
formation. □ 

Example 22. Consider the quantum graph shown in Figure [7} This graph has 
one vertex, one open edge, and one loop edge. Assume that the length of the loop 
is L. We examine the sheaf on three different open sets, U , V , and W as shown 
in the figure. The dimension of the space of sections over U and W is 2, and the 
dimension of the space of sections over V is 3. In the figure, basis elements (in 
terms of a traveling wave decomposition ) are shown. 

If we consider sections that are defined over UUV , they must agree on U n V. This 
induces the following gluing conditions: u\ — v\ and u-i — \v\ + |e lfcL W2 + |e ifcL i>3. 
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Observe that since V and W are disjoint, there are no gluing conditions required 
to construct sections over VUW . If instead, we moved W to W , so that VClW =/= 0, 
two gluing conditions apply: w\ — U3 and W2 = §i>i + ^e lkL V2 + ^e lkL Vs. 

4.3. Computation of the space of solutions. The objective of this section is 
to completely characterize the cohomology of transmission line sheaves in terms of 
the geometry and topology of the underlying metric graph. The main result is the 
following theorem: 

Theorem 23. Suppose Q is an F ' -transmission line sheaf with on a connected quan- 
tum graph that has I closed edges (of which V of them have edge endomorphisms E 
that satisfy E — E^ 1 = 0, a resonance condition for closed edges), m open edges, 
and n resonant loops. Then 



(9) dimH°(X;g) 




(10) dimH l {X-g)^ 



if I = m = 
i/! = 0,m/0 
min{0,Z'-l} ifl^0,m = 
min{0, /' — 1} otherwise 



n + 1 if I = m = 

if I = 0,m ^ 
+ 1 + min{0, /' - 1} ifl^0,m = Q 
+ min{0, 1' — 1} otherwise 



This rather long section begins first with a discussion of some of the implications 
of this result in Section l4.3.1l The proof of Theorem 1^31 is outlined in Section [4.3.21 
with the proofs of the key edge collapse lemmas being postponed until Section l4.3.3l 

4.3.1. Sheaf obstruction theory for quantum graphs. Of course, Theorem l23l stronglv 
restricts the geometry of the underlying graph, but it also restricts the topology. 
In this section, we therefore also answer the question: 

Question 24. Given the cohomology of a transmission line sheaf J 7 on a quantum 
graph X, how much does it restrict the topology of X? 

If there is loss in the quantum graph then no loop and no closed edges are res- 
onant. Therefore, in lossy quantum graphs, the cohomology is entirely determined 
by the number of open edges and whether there are any closed edges (resonant or 
not). The space of excitations of a lossy quantum graph is essentially determined 
by external scattering effects only. 

Clearly, a graph with no resonant loops and no resonant edges will have trivial 
H 1 . The converse is not true, though. Closed edges will typically induce a trivial 
H 1 if considered on their own, regardless of edge endomorphism. This means that 
the cohomology of F-transmission line sheaves is not a strong enough invariant to 
distinguish between graphs that contain closed edges and those that do not. In 
particular, consider a graph X that consists of one vertex and n open edges, and a 
graph Y that consists of one vertex, n open edges, and one closed edge. Regardless 
of the edge endomorphism for the closed edge in Y , any transmission line sheaf 
over X will have the same cohomology as any other transmission line sheaf over 
Y, though clearly their geometry differs. This is a sheaf-thcorctic expression of the 
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existence of isospectral quantum graphs j35], but goes a little farther: there exist 
isospectral graphs whose eigenspaces agree as well. 

A convenient way of summarizing the above discussion is that nontrivial sheaf 
cohomology classes constitute obstructions to the topology of the graph being triv- 
ial. 

Lemma 25. Suppose that J- is an ¥ -transmission line sheaf on a metric graph X 
which has no closed edges. Then 

• the number of open edges in X is equal to the Euler characteristic of JF , 
namely x(^ r ) = dim H®(X;J-) — dim H 1 (X;J-), and 

• there are at least dim J 7 ) loops in X. 

The appearance of the Euler characteristic in Lemma[25]fits neatly with the inter- 
pretation of H 1 as describing resonances. In particular, one can think of the number 
of open edges being the total number of excitations (dim H°(X; J 7 )), which includes 
both externally-induced excitations and internally-induced excitations (resonances) 
minus the number of resonances (dim H 1 (X;J 7 )). Therefore, ihe cohomology of 
the excitation sheaf over a metric graph encompasses both scattering problems and 
eigenvalue problems. 

If we consider the more specific case of lossless excitation sheaves, then the exact 
number of loops can be determined. Suppose X is a finite metric graph with no 
closed edges. As noted in Lemma [25j the number of open edges can be computed 
from the cohomology of any excitation sheaf. Each edge of the graph has a (set of) 
resonant frequencies, some of which may coincide. Generically, none of the resonant 
frequencies coincide, and it is therefore easy to locate the lowest resonant frequency 
of each edge, which determines the length of each edge. This is essentially the idea 
of [T7]. However, the genericity condition can be relaxed, since the cohomology of 
excitation sheaves captures edge multiplicity information. 

This discussion can be summarized by the following proposition and algorithm. 

Proposition 26. The topology and edge lengths of a finite metric graph with no 
closed edges is completely determined by the cohomologies of the excitation sheaves 
over it. 

Algorithm 27. Suppose that X is a finite metric graph with no closed edges, and 
the £k is the excitation sheaf on X with wavenumber k. 

(1) Determine the number of open edges m using Lemma Wb\ 

(2) Define 

_U{k,p)\keM+ and dimH 1 {X;£ k )=p^Q} ifm^O 
° ~ \{(k,p)\k e R+ and dim if 1 (X; £ fc ) = p - 1 + 0} if m = 

This is the set of resonant wavenumbers of loops in X, counted with multi- 
plicity. 

(3) Compute hi = min{k\there exists ap 6 N such that (k,p) £ Si}. This is 
a fundamental resonant wavenumber of an edge in X . Notice that this is 
well-defined since Si is countable. 

(4) Define 

Si+i = {(k,p) S Si\k £ hN} U {{k,p- l)|(fc,p) e Si and k e hN andp> 1} 

(5) Iterate steps 3 and 4- Notice that since X is finite, the Si stabilize at the 
empty set after finitely many iterations. 
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Figure 5. Sets that cover the edge e which is to be collapsed 



4.3.2. Proof of Theorem\2^ The proof relies on three edge collapse results (Lemma 
1291 Lemma l3"0I and Lemma [3~T1) that permit combinatorial simplification of the quan- 
tum graph without disrupting the structure of the solutions. Combinatorial edge 
collapse is not new, and plays an important role in the quantum graphs literature, 
for instance |26j . These lemmas permit a direct, explicit computation of the co- 
homology of transmission line sheaves. (For a similar computational methodology, 
which finds a natural incarnation as the sheaf theory we mention here, see 0.) 

Definition 28. Suppose that X and Y are finite metric graphs, which may be 
disconnected. A map f : X — > Y is called an edge collapse if 

(1) There exists an edge e of X for which /(e) is a vertex ofY, 

(2) / restricted to X — e is a homeomorphism onto its image, which is Y — f(e). 

The first edge collapse lemma applies to a flow sheaf on a directed graph, and is 
the most general result for these kind of sheaves. 

Lemma 29. Suppose that f : X —± Y is an edge collapse of an edge e that has 
distinct endpoints, and that T is an ¥-flow sheaf on X . Then f induces an isomor- 
phism on cohomology: 

(11) H*(X;T) = H*(Y;.f*T), 

and f*T is an F-flow sheaf on Y in which the coding maps at the vertices in Y— /(e) 
are unchanged from those in X — e. The coding map at e is given by the following. 
Let V\ and v 2 be the endpoints of e, in which e is incoming for V2- Without loss of 
generality, suppose that e is the first output of v\ and the first input to V2 ■ Then 
(12) 

(f>i- 1 (v 1 )(a 1 ,a 2 , ■■■a m ) if i < n - 1 

^i- n (v 2 ){4>i{vi){ai,a2, ...a m ),b 2 ,b 3 , ...b p ) otherwise 

where the input degree of v\ is m, the output degree of vi is n, and the input degree 
of v w is p. See Figure^ for a graphical representation of this situation. 



^i(/( e ))( a i) fl 2, —a m , h, b 3 , ...b p ) 



It is of course tempting to wonder if multiple edges between a pair of vertices 
can be collapsed together, something which would be a bit stronger than Lemma 
l29l The success of this depends delicately on the coding maps. In the case of 
transmission line sheaves, edge collapse works as desired. 
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I closed edges 




n loops 



m open edges 

Figure 6. A graph with one vertex, I closed edges, m open edges, 
and n loops 

Lemma 30. Suppose that f : X — > Y is an edge collapse for an edge e with distinct 
endpoints, each of which have degree greater than 1. Then for any transmission line 
sheaf Q with field coefficients, f induces an isomomrphism on sheaf cohomology and 
the direct image f*Q is also a transmission line sheaf. 

Of course, Lemma [30] specifically excludes the case of collapsing an edge loop, 
since such an edge does not have distinct endpoints. The final edge collapse lemma 
permits loops to be collapsed provided they are not resonant. 

Lemma 31. Suppose J- is an ¥ -transmission line sheaf over a graph X , and f : 
X — > Y is a map between graphs that collapses a nonresonant loop. Then f induces 
an isomorphism on cohomology. 

With these lemmas stated (we prove them in Section |4 . 3. 3[) . we now address the 
problem of proving Theorem 1231 

Proof, (of Theorem l23[) Begin by collapsing out all nonresonant loops using Lemma 
|3"T1 Then, obtain a minimal spanning tree for X, and collapse each edge in the tree 
using Lemma [30] Thus the sheaf cohomology of Q can be computed by computing 
the sheaf cohomology on a graph Y with a single vertex, I closed edges, m open 
edges, and n loops as shown in Figure [6] 

Choose a good cover for this space by selecting U as a contractible open set 
containing the single vertex of Y. To each loop j, associate two additional open 
sets Kj, Hj homeomorphic to intervals, which complete the cover as shown in 
Figure [7] 

Note then that the Cech cochain complex is 

(13) > ¥ m+2n F 4« ]p2i F 6n W 2l > q 

For convenience, let us consider the case where I = 0, and let d — m + 2n, which 
is the degree of the single vertex of Y. Organize the coboundary map 8 so that it 



INVERSE PROBLEMS IN GEOMETRIC GRAPHS USING INTERNAL MEASUREMENTS 23 




has the following block form (columns are ci, 



where 
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We can do row-reduction on each Bq X 6 block individually, and in such a way that 
this preserves all other entries not in this block (in particular, the nonzero rows 



2 

d 

2 
d 
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remain unchanged outside the block). We obtain a new block 



/ 1 





-1 








°\ 


2-d 
d 


2-d 
d 














2-d 
d 


2-d 
d 

















1 











-1 








1 





-1 





V o 








1 





-v 



-^6x6 



which has rank 4 if d = 2 and rank 5 otherwise. Now observe that for the nonzero 
entries of the second rows of each B' block. They are 

2-d 2 2 2 

d d d "' d 

2 2-d 2 2 



2-d 
d 
2 

d 

2 
,/ 



2-d 



In each column, there are n — 1 copies of 2 /d and one copy of , so each column 



sums to 



1. Now when n 



m-\-2n 



ie. to = 0, this implies that the 



coboundary matrix has an additional kernel element, and clearly at most one such 
element. On the other hannd, if n ^ d/2, clearly all of these second rows of each 
B' block are linearly independent. Hence 



rank 5 = { 5n 

5ro 



if d = 2 

if TO ^ 

1 if to = 



Now if I > 0, so there are closed edges present, the coboundary map must be 
augmented with 21 rows and columns. Let d — 2n + m + l. For the i-th closed edge, 
we add two new rows, which look like 



1 

2-d 



-El 
Ei 



..m copies... 
..to copies... 



.2n copies. 
.2n copies. 





2 > 
d 



in which the first 21 columns are added, the next m columns are the first to columns 
of the original S, and the remaining columns count off in pairs from the first two 
columns of each Bq x q block. More precisely, they correspond to the columns la- 



belled ei,/i, 
becomes 



•••5 &i i fi j • ■ • ; &l j fli Cl j ■ ••) C-m ; d\ , 



After one row operation, this 



2 

Ei — E- 



..m copies... 
..to copies... 



...2n copies. 
...2n copies. 



Note that we thereby obtain a duplicate copy of the second row for each closed 
edge, so at most 1 is contributed to the rank by these rows. One the other hand, 
the first row is clearly linearly independent from all the others. Hence the rank of 
the coboundary map is increased by max {21, 21 — I' + 1}. □ 

4.3.3. Proofs of the edge collapse lemmas. We begin by addressing the most general 
edge collapse result, Lemma The central difficulty is that the edge endomor- 
phisms and coding maps are not specified with a particular form. This complicates 
the calculations somewhat. 
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Proof, (of Lemma |29|) We aim to employ the Vietoris mapping theorem to obtain 
the desired isomorphism on cohomology. To this end, observe that since / is an 
edge collapse, it follows that it is a closed surjection. Additionally, X and Y are 
both paracompact, so / _1 is always taut. Suppose that y € Y, and discern two 
cases: 

(1) That y is not /(e), in which case / _1 is exactly one point, so H p (f~ 1 (y); F) = 
for p > 0. 

(2) If y = /(e), observe that H p (e; J 7 ) = lim H p (U a ; J 7 ), where U a ranges over 
open sets containing e. We consider a good cover of U a that consists of V\ 
(containing the vertex v\) and V2 (containing V2) whose intersection lies in 
the interior of e. The Cech complex is then 

( 14 ) > M m © M p — S —> M > 

In which the coboundary map 3 is given by [a%, 02, ■■■a m , b\, 62, ---bp) <— > 
4>i(vi)(ai, C12, ■■■a m ) — b\. Since 4>i( v i) is a homomorphism, it's clear that 
the image of 8 is F. Hence, H p (e; J-) = for p > 0. 
For the second statement, observe that the only thing to check is that the stalk 
over /(e) has the correct rank. In this case, that rank is m-\-p — 1, which agrees 
with Definition [H □ 

Remark 32. The formula for the coding map at the collapsed vertex /(e) can be 
written in terms of matrices (with entries in ¥). Suppose A is a matrix for 4>i v i) 
and B is a matrix for ${1)2) ■ Let u T be the first row of A, which corresponds to 
the output of vi along the edge e. Likewise, let v be the first column of B, which 
corresponds to the input to V2 coming from e. Let a and b be the matrices obtained 
by deleting the first row of A and first column of B, respectively. Then </>(/(e)) has 
block matrix form 



(15) 0(/(e)) 



a 
m T b 



Corollary 33. If J- is an F-flow sheaf over a finite, connected graph, its sheaf 
cohomology can be computed by looking at the direct image under the collapse of a 
spanning tree. 

For the proof of Lemma 1301 in addition to verifying that the Vietoris mapping 
theorem still holds, we must also verify that the direct image is still a transmission 
line sheaf. This requires a straightforward, but lengthy, computation. 

Proof, (of Lemma l30|) We need only redo the case of H* (/ _1 (/ (e)); J 7 ) in Lemma 
I2U1 In this case, look at a good cover of the edge e, which consists of three sets 
{U±, U2, U3}. Let U\ contain one endpoint of e (degree m), U3 contain the other 
(degree n), and Us lie entirely within the interior of e. We'll assume that the e has 
an edge endomorphism L. 

In the Cech cochain complex, the coboundary map has the form 

°\ 

-1 



which we claim has rank 2. (The columns are organized by G (Ui) ® G (U3) ® G (U2) .) 
Hence the H 1 ^^) = 0. 
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We now address the claim by examining the kernel of the coboundary map. By 
a pair of linear combinations of rows of (|16|) . we obtain 



\ m £ — ' 

and 



2 m_1 

Cj - c 

m ^ 

3=0 



2 11-1 

11 3=0 



Then we can solve for do by substitution: 



m 

j=i v 7 \ 3=0 

2 ^ 2/2 \ ^ / 2 \ / 2 



whence 



to J n \m J f-J \m J \n I 



m / \n J J m n \ m J 

_ . \ _ , 771—1 „ / „ \ n— 1 

2m + 2n - 4\ , 2L 2/2-m\v^, 

777 ^ 77 V m / ^— ' 



r m— 1 n — 1 

77_L v—v 2 — 777, 



( l7 ) rf o = ^ > c ? H > d,-. 

777 + 77 - 2 <^ 3 777 + 71 - 2 ^ J 



In particular, this confirms that the rank of (|16[) is 2, since the kernel is of dimension 

777 + 77 — 2. 

Continuing with (|17[) , we show that the direct image is a transmission line sheaf 
by exhibiting a typical restriction map to an output of one of the edges. Without 
loss of generality, we consider the output along the i-th edge of V2, namely 



2" ' 



n 

3=0 



d'i = - d 3 ~ d * 

E c * + i. ( 9 ) S + ^ S - * 



2L ^ 2 / 2 - 777 \ ^ . 2 " 1 



777 + 77 — 2 77 V 777 + 77 — 2 / 77 

i=i x 7 i=i i=i 



' m — 1 n—1 , 

E / - r - • E ) ~ di 



777 + 77 — 2 

which is of the form required for a transmission line sheaf. □ 



INVERSE PROBLEMS IN GEOMETRIC GRAPHS USING INTERNAL MEASUREMENTS 27 





LoL\ 



edge collapse 




ioL 3 



Figure 8. Edge collapse results in edge endomorphisms being composed 



Remark 34. Notice that when the edge e is collapsed using Lemma [WA the edge 
endomorphisms on one side are composed (see Figure^). This composition is non- 
unique: in the proof, we composed L with the edges on the side of V\, but we could 
have composed by L^ 1 on the side of v^. This clearly results in quasi-isomorphic 
sheaves, as they will agree in cohomology, and this difference will not concern our 
discussion here. 



Remark 35. We note that a degree 1 vertex causes the proof of Lemma [5U\ to fail 
because the direct image is generally not a transmission line sheaf. The Vietoris 
Mapping theorem applies perfectly well in this case, but the coding maps do in fact 
change. It is true that H 1 (e) is still trivial, even if there is a nonidentity edge 
endomorphism. 

Suppose for instance that m = 1, and that we wish to compute the restriction to 
the output of the i-th edge incident to V2- This has the value 



d' 



-£<*; 



If it the collapse of the 
have 



resulted in a transmission line sheaf, then we should 



d' = 



n-1 



1 ^ 



upon eliminating do using the conditions at v\. However, what we instead obtain is 

2 / L + L- 1 



d' 



(2/n- l)L + L~ 



n-l 



where L is the edge endomorphism for e. Although certain values for L will result 
in a transmission line sheaf, generic values of L will not. 

An interpretation of this result is that the degree 1 vertex's influence is to adjust 
the coding maps, essentially "tuning" the transmission line. In the case where a 
transmission line sheaf results from an edge collapse of a closed edge, the edge is 
the correct length to have no effect at all, which is related to resonance phenomena. 
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a 



degree k vertex 




Figure 9. Cover of a loop in the graph, for use with the Mayer- 
Vietoris sequence 

Finally, we address the case of Lemma ED Rather than using the Cech ap- 
proach as in the earlier calculations, we instead use a Mayer- Vietoris sequence, to 
illustrate an alternate technique for cohomology computation. This has the advan- 
tage of requiring fewer dimensions, and is a fairly natural context to consider the 
computation of sheaf cohomology over graphs. 

Proof, (of Lemma UTT]) Let U be a connected open set covering the loop to be 
collapsed, and contain exactly one vertex of degree k. We form a cover of U by two 
open sets A and B. Let A be a connected open set contained in the interior of the 
edge, and B be a connected open set containing the unique vertex in U (see Figure 
U). Observe that IF (A) = F 2 and F(B) = ¥ k . The appropriate Mayer- Vietoris 
sequence is 



H°{U;T) ^ H°{A:T)®H°{B;T) — - — ► H°(An B; J 7 ) — - — y H 1 (U;J : ) ->■ 0, 

where we observe that H°(A R B; IF) = F 4 . As a result of this exact sequence, 
H a (U;F) = ker A, and H 1 (U;J') = image S. Observe that 



Suppose that the edge endomorphism for the edge to be collapsed is L : F — >• F. 



dim ker S = 



4 - (k + 2) + dim ker A 
2 - fc + dim H°(U; T). 



Then 




using Kirchoff conditions at the vertex. Performing some algebraic manipulations, 
we find that elements of the kernel of A satisfy (provided k ^ 2) 
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and 



= ( I - \lA b +fi-^jLa-^(c 3 + ... + c k ) . 



If k ^ 2, then this leads to 

= (2 - L - L- V )b +(L- l)(c 3 + ... + c fe ), 

which means that the dimension of the kernel of A is k — 2 if L =/= 1 and k — 1 if 
L = 1. 

If instead, k = 2, then we obtain 

a = La, = L~ l b, 

which implies that the dimension of the kernel of A is if L ^ 1 and 2 otherwise. 
Thus we have 




if L = 1 and k ^ 2 
if L ^ 1 and k ^ 2 
if L = 1 and fc = 2 
if L ^ 1 and k = 2 



and 



'F if L = 1 and fc ^ 2 

if L ^ 1 and k ^ 2 

F if L = 1 and k = 2 

^0 if L ^ 1 and fc = 2. 

The conclusion is that the Vietoris Mapping theorem applies for the loop to be 
collapsed if and only if it is nonresonant, that is, if L ^ 1. □ 

4.4. Geometry extraction algorithm. If the topology is known, then generically 
a single section of the excitation sheaf contains all of the geometric information. 
This result is stronger than the results obtained by previous authors (like Proposi- 
tion [26]) which requires knowledge of many sections. However, the sheaf-theoretic 
framework provides a local, iterative mechanism for describing the geometric in- 
formation in a quantum graph. The central idea is that in collapsing a spanning 
tree in the graph, the cohomology (and hence the sections) of flow sheaves is un- 
changed. Under such a collapse, however, the global sections are very easily tied 
to the metric structure of the graph. By sequentially "undoing" the edge collapses, 
edge endomorphisms are determined one at a time until all are determined. 

Operationally, enough information can be obtained by placing a directional sen- 
sor at each vertex with degree not equal to 2. Each such sensor detects the incoming 
wave amplitude along each incident edge. Since the algorithm measures phase dif- 
ferences between points in the graph, this requires that the sensors be synchronized, 
and take their measurements simultaneously. From a practical point of view, loss 
in the graph (which we have neglected in this section) limits the visibility of the 
signal sources. As such, it is probably unecessary to require synchronization over 
all sensors placed in a lossy graph. 

Theorem 36. Suppose that the edge endomorphisms of a transmission line sheaf 
T are algebraically independent (in the ring of ¥ -endomorphisms). In the case of 
quantum graphs, this is equivalent to requiring that the edge lengths are algebraically 
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independent. Then the edge endomorphisms are determined by any nonzero section 
of J 7 when Kirchoff conditions ((8J) are used. 

As might be expected, the proof of this result is more interesting than its state- 
ment, and proceeds by an inductive computation. 

Proof. Base case: We assume that X consists of a boquet of circles, namely that 
there exists a single vertex v in X; we determine all of the edge auomorphisms in T 
from s. For concreteness, assume that there are n closed loops and m open edges in 
X. Consider a connected open set U which contains the unique vertex and none of 
loops completely. Then T(U) will be isomorphic to F 2n+m , that is there are 2n + m 
incoming signals entering the vertex. Some of these, of course, will be related if 
we consider T(X). Suppose without loss of generality, that dx, g?2, d^n-i, are 
the values of the section s on the loops at v, and that ci, ...,c m are the values of 
the section s on each open edge at v. We can further organize the loops so that dj 
and dj+i are the values at either end of a loop, following Figure |5J Hence, we can 
use ((SJ to obtain a set of n equations, one for each loop 

which can easily be solved for the edge endomorphisms Lj. 

Inductive step: We assume that there is a graph Y which can be obtained 
from X by collapsing f : X ~^Y a single edge e with distinct cndpoints v%, v%, and 
that all edge endomorphisms of f^J 7 are known. Assume (see Figure [5]): 

• That the degree of v\ is m and that the degree of v<i and n, 

• That the edge endomorphism of e is L, 

• That Co, c\, c TO _i are the values of the section s, that are incoming for 

• That Co is the value of s coming into V\ along e, and 

• That do, di, d n _i are the values of s incoming to 1)2, 

then ([8]) gives the following equations (see Figure [5] for notation) : 

2 2 
(18) b = — (c + ... + Cm-i) - c , a = ~(d + ... + d n -i) - ^o- 

m n 

Using the edge endomorphism i, we note that Lb = do and Co = L~ 1 a, which 
permits each equation in (|18[) to be solved for do or L. 

d = — (c + ... + c TO _i) - Lc = c - (di + ... + d„_i). 

m 2 — n I — n 

Since the graph X is finite, it only remains to see that the induction can be 

started by finding a sequence of trees = T± C T2 C ... C T in which T, — Tj_i is a 

single edge, and T is a spanning tree for X. The base case is obtained by examining 

/ : X — > X/T, and each induction step is obtained by considering the collapse via 

/, : X/Ti-i X/Ti. In each case, these maps satisfy the hypotheses of Lemma l30l 

□ 

5. Discussion 

In order to implement our algorithms into a viable sensing system, one needs 
to address discretization issues. For instance, the topology extraction algorithm 
assumes that there are enough receivers to discriminate whether a coverage region 
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(or intersections thereof) is disconnected. This may require a very high density of 
receivers in order to make this discrimination with confidence, since it essentially 
amounts to measuring the distance between clusters in signal space. Although such 
discretization effects are out of the scope of this article, it is useful to note that the 
amount of loss, the distribution of edge lengths, and the operating frequency play 
an important role in determining the necessary receiver density. 

As a related point, detecting the appropriate threshold to use for the visibility 
regions is less clear when using discrete receivers. A rigorous approach to this 
problem might use the tools of persistent homology [TU] to attempt to capture the 
topology of the visibility regions, and select appropriate thresholds. 

While the geometry extraction algorithm assumes discrete receivers, their syn- 
chronization constitutes a major limitation to performance. Indeed, the synchro- 
nization requirement is modulated by the mutual visibility of the receivers, which 
depends on signal loss. On the other hand, the synchronization requirement may 
also be relaxed by the use of active sensing, in which the signal sources and sensors 
are colocated. This method would suggest that instead of studying homogeneous 
solutions, we instead study fundamental solutions. 

Finally, in practical urban imaging applications, the graph model is inaccurate 
for open areas, though in simulation sensible results can be obtained [T3]. Indeed, 
one needs to generalize the theory discussed here to handle higher-dimensional cel- 
lular spaces. Many difficulties arise from this generalization, not the least of which 
is that the cohomology of the excitation sheaf becomes infinite-dimensional. The 
refinement algorithm as stated here no longer works correctly, either, as connect- 
edness no longer implies contractibility of intersections. Thresholding to obtain 
contractibility should still work correctly, though verifying that it does will require 
considerable effort. 
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